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Anomalies near the conductance threshold of nearly perfect semiconductor quantum wires are 
explained in terms of singlet and triplet resonances of conduction electrons with a single weakly- 
bound electron in the wire. This is shown to be a universal effect for a wide range of situations 
in which the effective single-electron confinement is weak. The robustness of this generic behavior 
is investigated numerically for a wide range of shapes and sizes of cylindrical wires with a bulge. 
The dependence on gate voltage, source-drain voltage and magnetic field is discussed within the 
framework of an extended Hubbard model. This model is mapped onto an extended Anderson 
model, which in the limit of low temperatures is expected to lead to Kondo resonance physics and 
pronounced many-body effects. 
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I. INTRODUCTION 

Semiconductor quantum wires can be fabricated with 
effective wire widths down to a few nanometers; for ex-, 
ample, byJieteroepitaxial growth on iv'-groove surfacesu 
and ridgesa, cleaved edge over-growthp, etched wires with 
gatingp, and gated two-dimensional electron gas (2DEG) 
structureg3'El. More recently, there has been consider- 
able interest in carbon nanotubes for which the quantum 
wire cross section can approach atomic dimensions. Such 
structures have potential for opto-electronic applications, 
such as light-emitting diodes, low-threshold lasers, single- 
electron devices and quantum information processing. 

Conductance steps in various types of quantum point 
contacts and. quantum wires were found more than a 
decade agcOu. These first experiments are broadly con- 
sistent with a simple non-interacting pictureQ. However, 
there are certain anomalies, some of which are believed 
to be related to electron-electron interactions and appear 
to be spin-dependent. In particular, a structure is seen 
in the rising edge of the conductance curve, starting at 
around 0.7(2e 2 //i) and merging with-, the first conduc- 
tance plateau with increasing energyQ. ■-■ This structure, 
already visible in the early experimentsB, can survive to 
temperatures of a few degrees and also persists under 
increasing source-drain bias, even when the conductance 
plateau has disappeared. Under increasing in-plane mag- 
netic field, the structure moves down, eventually merging 
with the e 2 /h conductance plateau at very high fields and 
is not a transmission effect through a ballistic channolc. 
A structure is seen also in high quality quantum wiresli^ . 
In some experiments, an anomaly is s^ee at lower energy 
with conductance around 0. 3(2e 2 /hptl. This can also 



survive to a few degrees, though is less robust than the 
0.7 anomaly and is more readily suppressed by a mag- 
netic ficldo. Recently the anomaly-, was confirmed also 
in back-gatedta, in shallowj£.tchedll3 point contacts and 
in a ballistic quantum wiretj. At low temperatures the 
anomaly exhibits a puzzling similarity with Kondo reso- 
nance behaviorEj, as do thermopower measurementsEEI. 

Theoretical work has attempted to explain these ob- 
servations in various ways, including conductance sup- 
pression in a-Xuttinger liquid with repulsive interaction 
and disordeiO, local spin-polarized density-functional 
theoryliS and spin-polarized sub-bandslij. Near the con- 
duction threshold, there is a 'Coulomb blockade' and we 
have shown that this gives rise to spin-dependent resp^ 
nances, for wires of both rectangular! 2 !! and cylindricalcil 
cross-section, with related anomalies in thermoelectric 
transport coefficientsEHl. A similar singlet-triplet sce- 
nario was presented in Ref. ^3] and a phenomenologi- 
cal approach is presented in Ref. |24|. Recent studies 
have investigated the 0.7 anomaly in quantum, point 
contacts within the Hartree-Fock approximationt3, spin- 
fluctuation backscatteringc3 and in the framework of 
the Anderson model with related Kondo resonance 
behaviorEZl. 

In Refs. p0|j2l| , p2] we suggested that these anomalies are 
related to weakly bound states and resonant bound states 
within the wire. These would arise, for example, from a 
small fluctuation in thickness of the wire in some region 
giving rise to a weak bulge. If this bulge is very weak 
then only a single electron will be bound. We may thus 
regard this system as an 'open' quantum dot in which the 
bound electron inhibits the transport of conduction elec- 
trons via the Coulomb interaction. Near the conduction 
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threshold, there will be a Coulomb blockade and we show 
below that this also gives rise to a resonance, analogous 
to that which occurs in the single-electron transistors. 
This is a generic effect arising from an electron bound in 
some region of the wire and such binding may arise from 
a number of sources, which we do not consider explicitly. 
For example, in addition to a weak thickness fluctuation, 
a smooth variation in confining potential due to remote 
gates, contacts and depletion regions could contribute 
to electron confinement along the wire or gated 2DEG. 
A significant contribution to the single-electron confine- 
ment could also arise from its electronic polarization of 
the lattice or image charge. 

In this paper we extend our previous study of a partic- 
ular geometry of the quantum wire with a comprehensive 
analysis of a wide range of shapes and sizes of wire in or- 
der to demonstrate the generic and wide applicability of 
the phenomena. We study in particular the threshold 
of the conductivity of nearly perfect wires for which a 
single electron is bound. We express the conductance 
in terms of the two-electron scattering matrix. In order 
to extend the exact two-electron analysis into the true 
many-electron domain, we construct an extended Ander- 
son model and analyze the influence of the corresponding 
momentum dependent coupling matrix elements. 

The model is introduced in the next section and the 
special case of a cylindrical GaAs wire is derived in Ap- 
pendix A. In Section III a detailed analysis of the two- 
electron problem is presented in which one electron is 
weakly bound in the wire, giving rise to spin-dependent 
scattering of the other. Exact singlet and triplet scat- 
tering states are computed near the conductance thresh- 
old. In Section IV we then show how the solutions of the 
scattering problem may be used to determine conduc- 
tance by an extension of the Landauer-Biittiker formula. 
This gives excellent agreement with a number of experi- 
ments on different kinds of quantum wire. The effect of 
finite magnetic field on the anomalies is presented and 
it is shown how they are related to the spin-split steps 
in perfect quantum wires. In the last Section we also 
examine the dependence of the anomalies on asymmetry 
introduced by finite source-drain voltage and summarize. 
Additional appendices are devoted to technical details on 
the solution of the two-electron wave function in an ex- 
ternal potential, and the Hartree-Fock analysis. 



II. BASIC MODEL 

In previous work@@il, we have considered a straight 
quantum wire with a small fluctuation in thickness giving 
rise to a weak 'bulge'. The precise details of the bulge 
are largely unimportant for what follows, the main re- 
quirement being that the change in the width of the wire 
is sufficiently gradual that inter-channel mixing of the 
transverse modes is negligible and that only one electron 
may be bound in the bulge region. The latter is always 
the case for a weak symmetric bulge, which has at least 



one bound state which can only sustain one electron due 
to Coulomb repulsion. The problem reduces to electrons 
moving in an effective weak potential well if we confine 
ourselves (by choice of gate voltage) to the Fermi energies 
for which no more than one transverse mode is occupied, 
i.e. the conductance threshold and the first conductance 
step. A typical effective potential well for such a bulge is 
shown in Fig. |l|. Such a potential well may arise in other 
ways, such as an actual potential fluctuation due to a 
nearby unscreened charged impurity, or even some self- 
consistent effect due to the electrons themselves through 
electronic polarization and image charge in a remote gate. 
We shall not consider the possible cause of this weak po- 
tential further but emphasize that because it may arise in 
many ways, the weak potential well model is very general 
with widespread applicability. 

e(z) 



z 

FIG. 1: Effective one-dimensional well caused by thick- 
ness fluctuation, impurity charge, gate image charge, self- 
polarization due to single electron or some combination of 
these. 

Consider now the motion of electrons in the wire near 
the conductance threshold. A single electron will be 
bound in the potential well region and the remaining 
electrons will undergo scattering from the localized elec- 
tron via the Coulomb interaction as they propagate from 
source to drain. At sufficiently low Fermi energy, the elec- 
trons in the source contact will be totally reflected by the 
bound electron due to Coulomb repulsion and there will 
be no current from source to drain at T = 0. As the 
Fermi energy is raised, the energy of the electrons in the 
source contact will be sufficiently high for them to over- 
come the Coulomb repulsion of the bound electron and 
a current will flow. In calculating this current we will 
make the approximation that the electrons flowing from 
source to drain only interact with the bound electron via 
a screened Coulomb interaction. This is a reasonable 
approximation provided that the electron density is not 
too low in the region of interest, i.e. the rising edge to 
the first conductance plateau. More precisely, the mean 
density of electrons in the wire (number per unit length) 
should be at least of order the inverse effective Bohr ra- 
dius of the material. We return to this point again in 
the final section. Within this approximation, the many- 
electron problem is reduced to an effective two-electron 
problem in which one electron is bound and the other 
is a representative electron at the Fermi energy in the 
leads. We show below that by solving this two-electron 
problem exactly and summing over all electrons near the 
Fermi energy we may compute the conductance. 
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A. Extended Hubbard and Anderson model 

The Hamiltonian corresponding to interacting elec- 
trons in the wire with a small geometric or potential 
inhomogeneity and close to the threshold of conduction 
is, within the effective-mass approximation an extended 
Hubbard Hamiltonian on a finite-difference latticed, 

H = H la + i ^2 UijUiUj + ^2 Uun^mi. (2.1) 

Here H\ a is the single-particle Hamiltonian: 

Hia = -t ( c l+i* c *<? + 4x c *+i CT ) + E ^ n i<y, ( 2 - 2 ) 



where c| CT , Ci a are electron creation, annihilation opera- 
tors, riia = c\ a Ci a , and = J2a ni <?- Model parameters 
are hopping i, local potential at site i, e,-, and screened 
electron-electron interaction at sites i and j, Uy. This 
Hamiltonian is derived and justified in Appendix A. 

In order to study the many-electron problem, it is 
also convenient to express the Hubbard Hamiltonian, 
Eq. (2.1), in a basis which distinguishes bound and un- 
bound states explicitly. Single-electron solutions corre- 
sponding to the tight-binding Hamiltonian Eq. (^^), fol- 
low from the single-particle Schrodinger equation 



Hi \<p) = E 1 \<p) , 



(2.3) 



and (with omitted spin index a), \tp) = J^j </?.j' c j |0)- ^or 
large |j| the potential 6j is constant, therefore the so- 
lutions are asymptotically plane waves. We thus diag- 
onalize this single-electron part of the Hamiltonian us- 



ing the transformation c qa = J2j Cja&p where 0* = 
(j\q) r~j exp(iqj) asymptotically for unbound states, with 
eigenenergies e q . In this basis the Hamiltonian becomes, 



H = J2 e i n i 



- E 
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W(9i92<7394) = Y Ui - 
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t>f. 



(2.4) 



(2.5) 



HamiltonianS0, 

H = e k n k + e d n d + y^(Vfe^dg4 ff d g + h.c.) + (2.6) 

Y2, Jkk'Sd ■ Skk'- 
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Un d ^n di + Y Mkk'n d cl a c k 'a 



kk> 



kk' 



Here U = U(dddd) is the Hubbard repulsion, V k = 
U(dddk) is mixing term, M kk , = U(ddkk') - \U{dkk ' d) 
corresponds to scattering of electrons and the direct ex- 
change cou pling is J kk i = 2U(dkk'd). Spin operators 
in Eq. (2/7) are defined as Sd = \ ^ aa , d\a a a'd a ' and 

Skk' = llZo-o-/ % a ^aa'Ch'a'i where the components of 
er are the usual Pauli matrices. A similar model has 
been proposed recently in Ref. Although the Hamil- 
toniarw-Eq. (|]^), is similar to the usual Anderson Hamil- 
tonianE2l, we stress the important difference that the kd- 
hybridization term above arises solely from the Coulomb 
interaction, whereas in the usual Anderson case it comes 
primarily from one-electron interactions. These have 
been completely eliminated above by solving the one- 
electron problem exactly. The resulting hybridization 
term contains the factor n dai and hence disappears when 
the localized orbital is unoccupied. This reflects the fact 
that an effective double-barrier structure and resonant 
bound state occurs via Coulomb repulsion only because 
of the presence of a localized electron. 

To be specific, we consider in this paper a cylindrically 
symmetric quantum wir£_with symmetry axis z and lat- 
eral coordinates r and </£j. Such a geometry corresponds 
to narrow 'v'-groove z-dependent quantum wires investi- 
gated recently, e.g. in Ref. O. The diameter of the wire 
is a (z), with zero potential within the wire and constant 
Vb > outside, i.e., 



V(r,z) 



0, r < \a (z) 
Vb, r > (z) ' 



(2.7) 



For the wire width, two generic shapes are taken, shown 
in Fig. @ with 



a (z) 
a(z) 



a o ( 1 — £ sin 7r- 



\z\ < ai 
\z\ > ai 



I ao (l + £cos 2 f^l , \z\ < ai 



, (2.8a) 



ao, 



\z\ > en 



(2.8b) 



We further denote the lowest bound state with energy 



e q < by d a 



-qtri 



with rid = Ylcr dtda and, similarly, 



the scattering states with positive e q are distinguished by 
q —>■ k. There are two independent unbound states corre- 
sponding to each k and these are chosen to be plane waves 

, . h 2 k 2 

asymptotically, i.e. $q — > e 3 as j — > ±oo and e k — ^jr- 
Retaining only those Coulomb matrix elements which in- 
volve both localized and scattered electrons, omitting all 
terms which would give rise to states in which the local- 
ized state is unoccupied, we arrive at an Anderson-type 




FIG. 2: The geometry of th e 'open quant um d ot' for the 
parameterization (a) Eq. (2.8a) and (b) Eq. (2.8b). 
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The region of interest is around z = and for large 
\z\ > ai the diameter is constant a . Single particle solu- 
tions corresponding to this geometry as well as the deriva- 
tion and calculation of parameters of the corresponding 
Hubbard Hamiltonian are presented in Appendix A. 



in the previous case and therefore the problem is closer 
to the strong interaction limit, as is seen also in Fig. ^, 
dashed-dotted line, where the dip in the electron density 
signals the strong interaction regime. The Hartree-Fock 
approximation gives too large energies here, which are, 
however, qualitatively correct. 



III. TWO-ELECTRON SOLUTIONS 

A. Bound states 

In order to calculate conductance though the system 
we first solve the two interacting electron problem for the 
present g eom etry using the extended Hubbard Hamilto- 
nian Eq. (2.1). Solutions for bound states are determined 
by numerical diagonalization of the system of equations 
presented in Appendix B, Eq. (B4). In Fig. || is shown 
the result of the two body electron density as a function 
of z/a for various shapes of the bulge [Fig. g(b)]. A 
general tendency is that long/narrow bulges correspond 
to stronger interaction resulting in formation of a dou- 
ble peak in density, as known from other studies of one- 
dimensional quantum dotaM As long as the two peaks 
are not well separated, the approximate methods men- 
tioned below are excellent, becoming gradually less reli- 
able with increasing separation between the peaks. 
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FIG. 3: Two-electron density for various bulge parameters. 

In Fig. |I] we present typical examples of the energy 
of two bound (singlet) electrons (E < 0) where 7 is the 
electron-electron coupling strength, defined by replace- 
ment U — ► ryU. Exact results are represented by the solid 
line, with other lines representing results obtained with 
the Hartree-Fock approximation, derived in Appendix C. 
At E > the lines correspond to the position of the 
singlet resonance, calculated with different methods and 
discussed below. In Fig. ^ the bulge is longer and nar- 
rower, therefore both singlet and triplet bound states ex- 
ist for small 7, while for stronger coupling the triplet is 
first pushed into continuum and finally, for 7 ~ 0.7, both 
states become resonances. Here approximate solutions 
are less accurate, because the bulge is much larger than 
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FIG. 4: Position of bound states (E < 0) and resonances 
(E > 0) vs. coupling 7, calculated exactly and within 
the Hartr ee-Fock approximation. Wire shape corresponds to 
Eq. ( 2.8b ) with parameters: ao — ai = 10 nm, £ = 0.24, 
Vo = 0.4 eV and k = 50 nm. 
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FIG. 5: As in Fig. ^ but with parameters: ao = 10 nm, 
ai = 4a , £ = 0.06, Vo = 0.4 eV and k = 50 nm. 



B. Scattering states 

Here we consider the scattering of an asymptotically 
free electron on a bound electron within the bulge. Such 
a system may be regarded as an "open quantum dot" in 
which one electron is bound and inhibits the transport of 
conduction electrons via Coulomb repulsion. The prob- 
lem is analogous to treating the collision of an electron 
with a hydrogen atom as, e.g., described in Ref. 32 and 
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studied by J.R. Oppenheimer and N.F. MottS. We only 
consider here cases in which the energy of the scattered 
electron is smaller than the binding energy of the bound 
electron. This ensures that only elastic scattering is pos- 
sible. Asymptotically the two-body wave function is a 
properly symmetrized product of a single particle bound 
state, \ip), and scattered state, \x(E)). 

For two electrons, the antisymmetrized wavefunction 
can be written as a product of a spin-part and an orbital 
part. We write the orbital part as 



4>i. 



Ipij + (-1) Ipji 

V2 



(3.1) 



ensuring that this is symmetric for singlets (S = 0) and 
antisymmetric for triplets (S — 1) (see Appendix B). For 
some large N 3> 1, ipij takes the form 



i < -N 
i > N 



(3.2) 



Asymptotic solutions for the unbound electron are ob- 
tained from the single-electron Hamiltonian Eq. (2.2) 
with the potential 



J2 a 



jk \Vk\ 



(3.3) 



for large \ Here \(p) is the single-particle bound state 
in the potential ej . Solutions with forward and backward 
currents have the following asymptotic form (j — > oo) 



Xj 



e ±ikj 

e ±i(kj-r,{k)ln2kj) K 



K < OO 
OO 



(3.4) 



for finite and infinite screening length k, respectively (see 
Appendix A). With, no screening (k — oo) Xj are the 
Coulomb functionso. 

Numerically exact solutions are obtained by solving a 
set of linear equations for the (2iV + l) 2 variables ifiij and 
transmission and reflection amplitudes. 



IV. CONDUCTANCE 

A. Single-electron solutions 

From the solution of the scattering problem, the con- 
ductance at zero temperature is .calculated using the 
usual Landauer-Biittikcr formalism^, 



G = GqT(E), 



(4.1) 



where Go = 2e 2 /h, E is the Fermi energy (in this case 
E = Ei) and T(E) is the total tra nsmi ssion probability. 

For an open bulge of shape Eq. ( 2.8a ), Fig. ||(a), thcx<i 
are no bound states and only single-electron solutionsL3 
are relevant. In Fig. ^ we present G as a function of 
electron energy for wires with shape Eq. ( p. 8a ) and three 



different widths. The main effect is a change of energy 
scale, according to scaling rule Eq. (A7), and the magni- 
tude of the conductance at the resonance energy, Go- In 
Fig. the cond uctance through the bulge of the shape 
from Eq. (2.8b) is presented. In contrast to the previ- 
ous figure, a bound state can exist here, indicated by the 
dashed vertical line. Further lines (n = 1,2) indicate 
bound states of individual channels for the special case 
when channel mixing terms in Eq. (A4) are set to zero. 
Dips in the conductivity in the second plateau (G ~ 3Go) 
correspond to Fano resonances caused by interchanncl 
mixing terms. 
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FIG. 6: G for a wire with shape Eq. ( |2.8a| ) and Vo = 0.4 eV, 
£ = 0.8, a\ = ao and in particular (a) ao = 7 nm, (b) ao = 
10 nm, and (c) ao = 15 nm. 




FIG. 7: G for wire shape Eq. ( p.8b| ) and Vo = 0.4 eV, £ = 0.5, 
ao = oi = 10 nm. Vertical lines indicate positions of bound 
states for the lowest channels. 

In Fig. ||(a) we again show the result of Fig. [j] compar- 
ing it with the one-channel approximation. In this paper 
we are interested in the rising edge of the conductance 
at the threshold, but with Coulomb interactions between 
bound and scattered electron included. Near threshold 
the one-channel approximation is excellent and therefore 
in the following we neglect higher channels. In Fig. 0(b) 
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is presented the influence of discretization parameter A, 
as introduced in Appendix A, on the conductivity. The 
position of the bound state is not strongly dependent on 
A (inset in enlarged energy scale) and for A < ai/5 the 
results obtained on the lattice agree with the continuum 
calculation within a percent, which justifies the use of 
the discretized Hamiltonian, Eq. (pTl|). 
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FIG. 8: (a) Conductance from Fig. ^ in comparison with the 
result obtained in the one channel approximation, (b) Con- 
ductance calculated with different A for wire parameters as 
in (a). Bound states are presented in the inset with enlarged 
energy scale. 



B. Interacting electrons 



We may extend the formula Eq. (4.1) to the case de- 
scribed in the previous section in which one electron is 
bound in the wire and the remaining electrons are trans- 
mitted with energy-dependent probability. Let P„ be the 
probability that the bound electron has spin a. It follows 
directly that the conductance due to all spin-up electrons 
in the leads is given by the extended Landauer Biittiker 
formula: 



G T = -[p T r TT (^) + p l r n (B)] 



(4.2) 



where Tf f is the transmission probability when the bound 
electron is spin up, Tj-j is the transmission probability 
when the bound electron is spin down and E is the Fermi 
energy. We have a similar expression for spin-down elec- 
trons in the leads and hence the total conductance is 



G(E) = - [P T T TT (i?) + PtTuiE) + P T T iT (£) + P L T U (E)\ 

(4.3) 

The transition probabilities 7tt and 7^ are different 
since in the former case the conduction and bound elec- 



trons both have the same spin (up) before and after scat- 
tering whereas in the latter case there are two possible 
final states, with or without spin-flip, i.e. 

T u (E) = \t n ^ n \ 2 + \t u ^ n \\ (4.4) 

where the scattering amplitudes are defined by 



<* T j I IV>u) 
(U 3 T \i>n) 



tu->nXi<Pj 



(4.5) 



as i — ► oo. IV'ti) is t ne exact scattering wavefunction 
and \iaja') — cj^cj^, |0). ipj — (J\(f) is the bound-state 
one-electron wavefunction and Xj = (j\x) is a forward 
propagating one-electron wavefunction at large j, as dis- 
cussed in Section III. 

In zero magnetic field it is clear that P-t = Pi = i in 



Eq. (4.3). We can express G(E) in a simpler form since 
the T aa i(E) are not all independent. Transforming to 
singlet and triplet base states (with S z = 0), 



\t,i,j) = 



NUI)-|H,jt) 

i»t,ji) + m,jt) 

V2 



(4.6) 



we get 



(Miil^u) 



t (0) 



(4.7) 



V2 



xm, 



where 

/"" - t 

fc u->T, 

Hence, Eq. ( |4~2] ) becomesEl 
1 



G — Gq 



(0) 



(4.8) 



(4.9) 



C. Results of numerical analysis 

In Fig. |s| we present the result of a comprehensive 
study of conductance for a variety of shapes of bulge. 
In Fig. ||(a) the bulge is wide and so short that only a 
singlet resonance is developed. The conductance there- 
fore exhibits a-structure similar to the 0.3 anomaly found 
in experiment^. In Figs. ^|(b,c) both singlet and triplet 
resonances are visible with a tendency for the resonances 
to sharpen as the bulge becomes weaker (£—►[)). 

In Fig. [h] the wire width is fixed at lOnm of the wire 
and positions of singlet (full lines) and triplet (dashed 
lines) resonances (or the corresponding bound states for 
E < 0) are plotted for various lengths and widths of 
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FIG. 9: Conductance for different shapes of the bulge. Each 
line is labeled with the parameter oi/ao. Other wire param- 
eters: ao = 10 nm, Vo = 0.4 eV and k = 50 nm. 



bulge, represented by ai/ao and £. We see that the reso- 
nances survive for a wide range of parameters. In Fig. |ll| 
is shown the position of singlet and triplet resonance en- 
ergies vs. the width of the wire, ao, with fixed shape 
of the bulge. The insets show the energy dependence of 
singlet and triplet transmission probabilities for selected 
special cases. Note that we have scaled the energy by 
a factor clqE. This would pro duc e identical curves for 
non-interacting electrons [Eq. (A7)]. 

After performing calculations for a wide range of pa- 
rameters, we conclude that a singlet resonance is always 
lower in energy than its corresponding triplet, in ac- 
cordance with Lieb-Mattis theorem- which however, is 
strictly valid only for ground states^. 
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FIG. 10: Energies of singlet (full lines) and tripl et (d ashed 
lines) resonances and bound states for wire Eq. (2.8b) as a 



function of ai/oo for different £. Other parameters are as in 
Fig. []. Full circles represent the energy, where the resonance 
energy is above the 'ionization' energy and the underlying 
wave function form is not valid anymore. 
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FIG. 11: The position of singlet (full line) and triplet (dashed 
line) resonances as a function of the width of the wire, ao. 
Note that the energy of the resonances is presented in a scaled 
form, ai = 50 nm, £ = 0.11 and other parameters are as in 
Fig- | 



D. Magnetic field along the symmetry axis 

The nature of conductance anomalies studied here can 
be further illuminated with experiments done in a strong 
magnetic fieldEJ. The effect of the magnetic field is, in 
our treatment, taken into account via the usual Zeeman 
splitting of channel energies. The incoming electron with 
Fermi energy E and spin component S z = ±| then has 
kinetic energy 



EtE b 



(4.10) 



where Eq = ^g* /j,bB. g* is the effective gyromagnetic 
ratio and /is is Bohr magneton. 



Near the conductance threshold we assume that the 
current is sufficiently low that the localized electron is in 
its ground-state with spin J, before each scattering event 
with a conduction electron. Hence Pi = 1 and P-t = in 



Eq. (4.2) which becomes 



GJE,B) = -T U (E,B) 



G l (E,B) = -T u (E,B). 



Since T n (E,B) 



(4.11) 



In this case Gi ^ Gt but, rather, 

(4.12) 

T n (£_,0) and T U (E,B) = 
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Tn(E + , 0), then the conductance is 
G = G| + Gj 

= UT n {E_,0)+T u (E + ,0)} 



(4.13) 



= e -[% (E+ , 0) + \% (E- , 0) + \% (E- , 0)] , 

where T s and T t are the same functions as in zero mag- 
netic field. 

In Fig. |l^(a) are plotted individual transmission prob- 
abilities for different spin configurations. Note that the 
spin-flip term is in general dominant at higher energies. 
In Fig. |lj(b) the corresponding results for conductance 
in the presence of a magnetic field is shown. The full line 
corresponds to B = 0, and other curves to B in incre- 
ments AB = 10T. 



1 

0.8 
0.6 

(N 

0.4 

0.2 


1 

0.8 

O 0.4 
0.2 




(a) 



(b) 



-4 



-2 




E (meV) 



FIG. 12: (a) Transmission probabilities for relevant spin 
configurations, (b) Conductance for B = (full line) and 
other lines for B in increments AB — 10 T. Parameters of 
the wire: ao = 10 nm, ai = 2.5ao, £ = 1.11, Vo = 0.4 eV and 
Hi = 50 nm. 



Note that Mkk' at higher energies tends to a constant, 
while other parameters approach zero, ensuring the cor- 
rect behavior at high-energy with unit transmission. The 
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FIG. 13: fc-dependence of matrix elements of the extended 
Anderson model. The wire is parameterized with ao = 10 nm, 
£ = 0.24, ai/ao = 2, Vb = 0.4 eV, k = 50 nm and 7 = 1. (a) 
Mixing coupling Vh- The energy + U is indicated with an 
arrow, (b, c) Scattering couplings M k y and J^k' ■ L is the 
length of the wire, where the wave functions are normalized. 



E. Results for the Anderson model 

As shown in Section II, the Hubbard model studied 
abo ve ca n be mapped onto an extended Anderson model, 
Eq. ( |2.7D . Conductance through a quantum dot described 
by a standard Anderson model is basically described by 
a peak or several peaks and at higher energies the con- 
ductance approaches zercH In the case of an open quan- 
tum dot, studied here, at higher energies the conductance 
tends toward unity, as a consequence of additional cou- 
pling parameters in the extended model. Here we analyze 
these terms individually and show their relative impor- 
tance. 

The coupling parameters are momentum dependent 
and in Fig. [13] the couplings Vfe, Mkk' and Jkk 1 are shown. 



scattering solutions of the Hamiltonian Eq. (2.7) are then 
obtained exactly for two electrons with the boundary con- 
dition that for z — > oo, one electron occupies the lowest 
bound state, whilst the other is in a forward propagat- 



ing plane wave state, <fik(z) 



^tkz 



From these solutions 



we compute the conductance again using the Landauer- 
Biittikcr formula. 

In Figs. |IJ(a,b,c) we compare the results of T s , T t and 
conductance G for a wire with the bulge as in Fig. |^. 
The thin lines are the exact scattering result for two elec- 
trons. The solid lines show the exact scattering solutions 
for the Anderson-type Hamiltonian, for which the matrix 
elements, and their energy dependence are calculated ex- 
plicitly. The solution of this Anderson-type model for 
two electrons, in which the localized level always con- 
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tains at least one electron, reproduce the main features 
of the exact scattering solutions of the original model. 
The energy dependence of the matrix elements is essen- 
tial to get this good agreement. Figs. |lj(d,e,f) show the 
corresponding results for a longer bulge from Fig. ||. Also 
shown in Fig. [li], dashed lines, are results with the direct 
exchange term omitted in Eq. (2.7). This term can have a 
significant quantitative effect, but does not qualitatively 
change the conductance curves. 

We have also solved a similar model in which plane 
waves, rather than exact scattering states of the non- 
interacting problem, were used. However, this gave poor 
agreement with the exact results. We conclude that an 
Anderson-type model is adequate for a near-perfect quan- 
tum wire provided that a suitable basis set is used and 
the energy-dependence of the matrix elements is accu- 
rately determined. Future work will focus on the many- 
electron properties of this effective Hamiltonian, includ- 
ing 'Kondo' and 'mixed valence' regimes. 




10 15 20 5 10 15 20 

E (meV) E (meV) 

FIG. 14: Singlet (a, d) and triplet (b, e) transmission prob- 
abilities and corresponding conductances (c, f). Parameters 
for the left set are as in Fig. for the right set: ao = 10 nm, 
£ = 0.15, ai/oo = 4, Vo = 0.4 eV, k = 100 nm and 7 = 0.9. 
Thin lines represent e xac t results from Eq. (2.1), thick lines 
are results from Eq. ( |2.7[ ). D ashed lines show results where 
the exchange term in Eq. (2.7) is neglected. 



F. Approximate methods 

It is not easy to get sufficiently accurate numerical so- 
lutions for the case of more than two electrons. Therefore 
it would be extremely useful if an accurate approximative 
method could be applied. The simplest approximation 
(presented here for the case of two electrons) is be the 
first iteration in solving the Hartree-Fock equations, as 
we presented in Section III for the case of bound states. 

We assume the two-body wave functions consist of a 
single particle state |y>^) and scattering state |x) with 
energy E. The two-electron wavefunction has then the 
form (see Appendix C), 



Ew*4 5A)t i°>- 



(4.14) 



The coefficients (fi are known, therefore only coefficients 
Xi must be determined. For the singlet state the simplest 
approximation is obtained if we perform the first iteration 
of the Hartree-Fock method subject to additional condi- 
tion that the electron has energy E ( "restricted Hartree- 
Fock" approximation): 



(Q\c i H 1 \x)+Y, U nWA 2 Xi=Exi 



(4.15) 



where, using Eq. (gj), (0 \ciHx\x) = -t{Xi-i + Xi+i) + 
CiXi- This is just the tight binding results for a single 
electron moving in an effective potential ej+X^j ^ij Wj : \ 2 ■ 
For the triplet state, the result is: 



(0\c i H 1 \x)+Y, U *i\'Pi\ 



Y^UtjtpjtpiXj = E X j- 
j 

(4.16) 



A better approximation for the singlet case starts from 
the unrestricted Hartree-Fock approximation, where the 
energy is 

{ip\H\i>) = (vj|JTi| <p) ( X |X) + <X X) (<p\<P) + 
+ (-if [{<p \Hr\x) (x\<p) + (x \Ht\<p) {<p\ X )) + 



^Uij \<piXj +xm\ 



2 ^ 

and the norm is given with 



«#} = £ \xi 



(4.17) 



Em* 



(4.18) 



The coefficients Xi are calculated from the Hartree-Fock 
equations based on the variation principle ( "unrestricted 
Hartree-Fock" approximation), 

(0 \c i H 1 1 x) + U v I^J I 2 Xi + U iiV*jViXj + 
j 3 
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In Fig. |l5| the effective one dimensional potential in 
Eq. fl4.15p is plotted for 7 = 0,0.5 and 1. The shaded 
region represents the position and the width of single 
particle resonance in this effective potential. This reso- 
nance corresponds to the singlet resonance presented in 
Fig. Jl6|(a) (dashed line), for wire parameters given in 
Fig. pjtogether with exact result, full line, and calculated 
from Eq. J4.15| ). 

We also show in Fig. |l^(a) the exact result and the 
corresponding result for an unrestricted Hartree-Fock 
scheme for which the wavefunctions of up and down spin 
electrons are different. As expected, the unrestricted 
method gives a more accurate result though both meth- 
ods reproduce the main features, the main discrepancy 
being an overall energy shift. Similarly in Fig. |l6|(b) we 
present the corresponding triplet resonance curve for pa- 
rameters from Fig. |[ Again, the overall agreement with 
the exact result is good apart from an overall energy shift. 




FIG. 15: Effective potential from Eq. (4.15) and for the wire 
with parameters as in Fig. H. 
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FIG. 16: (a) Singlet resonance (paramet ers fr om Fi g. ^[). Ex- 
act result and approximations of Eqns. (1.15) and ( 1.19 ) are 
shown, (b) Triplet resonance (para mete rs from Fig. Bh. Exact 
result and approximation of Eqn. (fl.lq) are shown. 



V. SUMMARY AND CONCLUSIONS 

We have shown that quantum wires with weak longi- 
tudinal confinement, or open quantum dots, can give rise 
to spin-dependent, Coulomb blockade resonances when a 
single electron is bound in the confined region. This is 
a universal effect in one-dimensional systems with very 



weak longitudinal confinement. The emergence of a spe- 
cific structure at G(E) ~ | anc j G ~ | =£- j s a conse- 
quence of the singlet and triplet nature of the resonances 
and the probability ratio 1:3 for singlet and triplet scat- 
tering and as such is a universal effect. A comprehen- 
sive numerical investigation of open quantum dots us- 
ing a wide range of parameters shows that singlet res- 
onances are always at lower energies than the triplets, 
in accordance with the corresponding theorem for bound 
statescS. With increasing in-plane magnetic field, the 
resonances shift their position and eventually merge in 
the conductance plateau at G ~ e 2 /h. With increasing 
source-drain bias we have shown why the higher triplet 
resonance weakens at the expense of the singlet, with the 
latter surviving to the point where the conductance steps 
themselves disappear. 

The existence of the conductance anomalies is a di- 
rect consequence of an effective double-barrier potential 
seen by the conduction electrons propagating from source 
to drain contacts under the influence of a bound elec- 
tron. For a symmetric one-electron confining potential, 
the existence of a bound state is guaranteed but this is 
not necessarily the case when the confinement is asym- 
metric. Such asymmetry in the confining potential may 
be easily achieved under a finite source drain bias and 
indeed, this was reportejdrip some of the experiments 
on gated quantum wiresocZI. These experiments show 
that as the source-drain bias is increased from zero, an 
anomaly appears at G ~ 0.25(2e 2 //i), coexisting with 
the 0.7(2e 2 //i) anomaly. Eventually, at larger bias, the 
remaining anomaly also disappears but only when the 
conductance steps themselves are on the point of disap- 
pearing, showing that the singlet anomaly is extremely 
robust. This behaviour is consistent with our model since 
under bias the triplet resonant bound-state will even- 
tually disappear because the confining potential in the 
x-direction will only accommodate a single one-electron 
bound state, giving rise to a singlet resonance only. This 
is shown schematically in Fig. |l7| where we also indicate 
the surviving singlet becoming broader with increasing 
bias resulting in a more pronounced step, as observed. 

Finally, we speculate on the exciting possibility that 
these anomalies in conduction are themselves a signa- 
ture for a new kind of conducting state in ultra clean 
wires close to the conduction threshold. Indeed, there is 
some experimental evidence for this in that the anomalies 
discussed above merge into a conductance step at e 2 /h 
under quite moderate magnetic fields and in the clean- 
est samples this behaviour is sometimes even seen in zero 
magnetic field. This suggests that there may be an under- 
lying spin polarized state associated with the propagating 
electrons in the quasi ID region. Such a spin-polarized, 
state would appear to violate the Lieb-Mattis theorerrc£l 
and would also need to be made consistent with our above 
explanation in terms of singlet and triplet resonances. In 
this respect we emphasize that the above theory must 
break down at very low electron density in the wire such 
that the mean separation between electrons in the wire is 
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FIG. 17: Effective double barrier showing singlet and triplet 
resonance with very small source-drain bias (a) and large 
source-drain bias (b). 



but not within the two-electron approximation we have 
used here and in some our earlier papers. It may well 
be that the two-electron approximation breaks down at 
low temperatures. The model presented here differs from 
the standard Anderson model in that the hybridisation 
term contains the factor n_ CT , and hence disappears when 
the localized orbital is unoccupied. This reflects the fact 
that an effective double-barrier structure and resonant 
bound state occurs via Coulomb repulsion only because 
of the presence of a localized electronj-JThe standard re- 
sults for the single impurity problem-^ thus cannot be 
applied directly to this effective model, and are a subject 
of current peascrchH However, a Kondo-like resonance 
is expectedE3'c3, for which many-body effects would dom- 
inate with a breakdown of our two electron approxima- 
tion. 



APPENDIX A: CYLINDRICAL WIRE 

1. Single electron basis 

A single electron in the wire considered here is de- 
scribed with the wave function ^(r, tp, z), which is a so- 
lution of the Schrodinger equation 



somewhat greater than the effective Bohr radius, the so- 
called strong correlation regime. In practical situations 
it is very difficult to avoid some kind of weak potential 
fluctuation which traps one electron. Indeed this may 
ultimately be impossible since even in a nominally per- 
fect wire, the presence of a single electron will polarise its 
environment leading to a potential well which will bind 
the electron giving rise to a Coulomb blockade for the 
remaining electrons, though the energy scale (tempera- 
ture) for this may be very low making it susceptible to 
masking the other effects. 

The main question is whether or not this confinement 
is sufficiently large for the electron density to exceed to 
the inverse Bohr radius when the wire begins to conduct. 
If the density remains low at this conductance threshold 
then we cannot ignore the mutual interaction between 
all electrons in the wire region, or even treat them self- 
consistently. In this situation, a more appropriate pic- 
ture would be one in which the Coulomb repulsion dom- 
inates and maintains roughly equal separation between 
the electrons as in a Wigner chain. On the other hand, 
if the mean electron density is of order, or greater, than 
the inverse Bohr radius, then an open quantum dot pic- 
ture with effective resonance levels for the propagating 
electron is more appropriate, as discussed in this paper. 

At low temperatures strong many-body effects are 
indicated frjQp, the activation-like behaviour jaf the 
conductanceBli3 and the thermopower coeficientO. As 
discused in our recend thermopower analysis, Ref. ^2], the 
anomaly at low-temperatures may well be a many-body 
Kondo-like effect contained within our extended Ander- 
son model, Eq. (pT7|), and studied recently in Ref. |27], 



2m 



■V 2 * (r, tp, z) + V (r, z) * (r, tp, z) = (r, tp, z) , 

(Al) 



where the effects of nonparabolicity are neglected and the 
effective mass is taken constant, m* = 0.067»*cioc with 
dielectric constant 12.5, appropriate for GaAsE3. 

At fixed z the wave function ^(r, tp, z) is expanded in a 
two-dimensional basis $ mn (r, tp; z) for the corresponding 
potential V (r; z) , Eq. (2.7). The coefficients in such an 
expansion over channels are ip mn (z) 

oc n 

#(r,<p,z) = ^ ^2 iPmn(z)$> m „(r,Lp;z). (A2) 

n— m— — n 

The transverse wavefunftions, $ mll (r, tp; z), depend only 
parametrically on z and take the form: 



(r, tp; z) = 



(z) Bin' (Kmn [z) r) + 
+Cmn (z) B m (n mn i z ) r ) 



GO 



'2m*e mn (z) 



'2m*\e mn (z)-V \ 



r < 
r > 



(A3a) 
(A3b) 

(A3c) 



°( z ) 

2 

a ( z ) 
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Here B, 



(i) 



7™. and B 



(2) 



eigenf unction: 
B 



K m are appropriate Bessel 
J m and 



(2) 



Y m , for e 



for e mn < V with B 

mn > U) ■ 



(1) 



The coefBcients A mn , B mn , 



Crnn and energies c mn are determined from the boundary 
conditions and the normalization of wave functions. 



Substituting Eq. (A2) into Eq. (Al) and integrating 
over r and tp leads to following coupled ordinary differ- 
ential equations for ip mn , 



\k 2 - k 2 

^ ^ bmnn 
n^n' 



(z) + a mnn (z)] ip 

■run + 

{ z Wm: 



E 



(A4) 
= 0, 



where the coupling coefficients are 



a m nn' (z) = 2tt / $* mn (r, v?; z) 

Jo oz 



^ o ^rnn' 
07 z 



f R d 

(z) = 4tt / <5>* mn (r, (p; z) — $ m „. 
Jo 



<9^ 



(r, </?; 2:) rdr, 
(A5a) 

(r, </?; z) rdr. 

(A5b) 



The coefficients coupling channels with different m are 
zero due to the orthogonality of e lmip for different m. 

Note that the Schrodinger equation, Eq. (Al), is in- 
variant under the transformation 



r 

E,V 



Ar, 

A~ 2 E, A~ 2 V. 



(A6) 
(A7) 



Extended Hubbard Hamiltonian 



We consider here the case when the variation in wire 
width is smal l, resulting in small derivatives of the coeffi- 
cients in Eqs. A5a , A5b . We consider only electrons with 
energy below the second channel and hence Eq. (A4) re- 
duces to a single equation for motion in z direction, with 
the potential 



e (z) = £oo (a (z)) + a 2 (z) e 00 (a (z)) 



(A8) 



The first term is the energy of the first channel and the 
second is related to aooo(-z) from Eq. (|A4|). a'(z) is the 
derivative of the wire diameter with respect to z. The 
second term in Eq. (A8) is always positive since 



e(z) -> e(»-e(oo) 



(A10) 



In Fig. |l8|(a) eoo (a) and eoo(o) are presented as a function 
of wire diameter. Figs. |l8|(b) and (c) show the variation 
of one-dimensional potential e[z) along the wire. 
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FIG. 18: (a) Dependence of eoo and eoo in Eq. ( |Aq ) on wire 
diameter for Vb = 0.4 eV. (b) One dimensional potential 
Eq. (A.S) for the wire shape Eq. (2.8a) and various values of 
£. Dashed lines correspond to the con tribu tion eoo- (c) The 
same as in (b) but for wire shape Eq. (2.8b). 



The single-electron Hamiltonian in the single channel 
approximation then becomes 



Bi = 



h 2 d 2 



<z). 



2m* dz 2 

This is readily generalized to many electrons, 



h 2 d 2 



V> CT (z) dz + 



(All) 



(A12) 



2m* dz 2 

+ \nih° w $1 w u & z ') ^ ^ dzdz ' 



eoo(a)= ^y 



& 2 tt [ R {d*oo(r,<p) 



da 



rdr. 



(A9) 



The potential Eq. (A8) is constant for large |z| and set 
to zero for convenience, i.e. 



where ipt (z) creates an electron with spin a at coordinate 
z and 



U(zi,Zj) 



47ree d(z i , Zj) 



(A13) 
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with 



1 = /' dridrj l*oo(r l ;z J )| 2 |$oo(r J ;^)| 2 



d(zi,zj) 



(A14) 



The Hamiltonian is further discretized at points Zj = jA, 
new creation operators are defined as 

4 = VA^t ( % .) . (A15) 
For sufficiently small A the difference formula is justified, 



dz 2 



Vv 0) 



i>a {Zi-l) ~ 2^<r (Zi) + Ipa {Zi+l) 



2A 2 



(A16) 



and Eq. ([All ) becomes the discretized extended Hub- 
bard Hamiltonian, 

H = H la + - ^2 UijUiiij + ^2 Uijn^Uji, (A17) 
where H\ a is single-particle contribution for spin c, 

Hla = -*y] (cj+l CT C «r + c L C i+l<?) + ^2 ei " io '' ( A18 ) 



with nj CT = cJg.Cio-, = "io-> hoping parameter, 



f = 



Tr 



2m* A 2 ' 



(A19) 



and Ci = 2t + e (zi). The effective distance between elec- 
trons at zi and za is after integrating Eq. (A13) over 
angular variables, 

87r / r i dr l I rjdrj |$ o (n; Zi)V x 



\/ {zi ~ Zj) 2 + {n - rj) 2 



where K is the complete elliptic integral of the first kind. 
d {zi, Zj) can be decomposed into distance along the wire 
and effective distance in the lateral direction, A (Zi, Zj) oo, 
i.e. 



(A21) 



d{Zl ' Zj) ^(^-Zjf + lXi^Zj)^}' 



The distance A ( zj, Zj) is invariant under the transfor- 
mation Eq. (|A6|), and hence the potential, Eq. ( |A13| ), 
transforms as 



For convenience we also take into account possible screen- 
ing with screening length k, i.e. 



U, 



Uije~ 



(A23) 



Under the transformation Eq. (A22) the screening length 
should be multiplied by A. In Fig. [l^ the parameter A 
is plotted for some typical cases, showing its dependence 
on wire width. 
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FIG. 19: (a) Lateral distance vs. separation along the wire, 
(b) Electron - electron interaction as a function of longitudinal 
separation. In both cases is ao = 10 nm and Vo = 0.4 eV and 
a =const. (c) Lateral distance at fixed z vs. wire diameter at 
Vo = 0.4 eV. 



APPENDIX B: TWO-ELECTRON WAVE 
FUNCTIONS 



Wave function for the case of two electrons are ex- 

(S S )t 

pressed in terms of a set of operators c^- ' creating an 
electron pair at sites i and j with spin S and z-component 
S~, i.e. 



|0). 



(Bl) 



U — > A~ 1 U. 



(A22) 



The base states 
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c g' 0)t 1°) 


c iT c il c il c iT 


10), 


(B2a) 


V2 


C g' m 10) 


= 4A io) , 




(B2b) 


cg' 0)t 10) 




10), 


(B2c) 




cg- 1}t 10) 


= 4A i°) • 




(B2d) 


form a complete set. 









If \tp) is a solution of Schrodinger equation 

H |V) = E |V) , 



(B3) 



then the coefficients Vij solve the system of linear equa- 
tions 



-lj + ipi+ij + tpij-i + ipij+i 
= (ei + ej + Uij - E) ipij, 
where we use compact notation 



In the basis Eq 



1 



(B4) 



(B5) 



the number of electrons on site i is 

2 



(VKI^) = 2 E r u ' 



the current for sites i and i + 1 is 



<V |/' i+1 |V) 



Im^ V*jV 



(B6) 



(B7) 



the energy is 

(V|#|V) 



should exist single particle bound states of the system 
and Ei is the energy of the state a. The energy of 
two electron states is shifted and defined to be zero, if 
one electron is bound and the other at the bottom of the 
single electron band, i.e. 



E^E-E- 



(i) 



(BIO) 



With this definition, the energy of two bound electrons 
is negative whereas it is positive when only one electron 
is bound. 



APPENDIX C: HARTREE-FOCK 
APPROXIMATION 

Here we neglect the Coulomb interaction between elec- 
trons (7 = 0). In the ground state both electrons are in 
the same state \ip) and the singlet wavefunction is 



IV) 



,(0,0)t 



V2 



|0) 



(CI) 



For finite 7 the best one-electron wavefunctions, (pi, are 
determined by minimizing the energy, 



cV* (VIV) 

Which leads to the equation 



= 0. 



(C2) 



d 



(V|ff|V)-(£ + £i 1} )^(VIV)=o, (C3) 



dip, 

where we have set 



(V|g|V) _ (1) 



(C4) 



taking into account the energy shift Eq. ( BIO ). The ex- 
pectation value of energy and the norm is 



> ' 7.1 



-^(e 4 + e, +(7„) 



Vi 



(B8) 



and 



and the norm of the wave function Eq. (Bl) is given with 

2 



(V|V) = {v\<pY 



(C5) 
(C6) 



(viv) = £ ft 



(B9) 



From Eq. (C3) follows a system of equations for coeffi- 
cients pi 



We consider quantum wires which are almost perfect 
but for which there is a very weak effective potential, giv- 
ing rise to bound states. The cross-sections of these wires 
are sufficiently small that the lowest transverse channel 
approximation is adequate for the energy range of in- 
terest. The smooth variation in cross-section also guar- 
antees that inter-channel mixing is negligible. We study 
here only wires with one weak bulge around z — 0. There 



(0\c l H 1 \p)+J2U iJ \p J \ 



EhfPi, 



(C7) 



where (0\ciHi\ip) — —t(ipi—i + p>i+i) + tip is a one- 
electron tight-binding Hamiltonian, y\ Tig \ipj\ is a 
Hartree potential and 



W -TP -I- (gjgjlg) 
— & + -C'l / 1 \ 



(C8) 
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is the so-called Hartree-Fock energy. The energy of a 
bound state is then given by 

E = 2E M - Uij \^\ 2 |^-| 2 ~ fif, (C9) 



where, due to double-counting of the interactions in sin- 
gle electron energies, Em is subtracted. In Fig. elec- 
tron density for singlet states with different 7 is presented 
for a shorter bulge with parameters as in Fig. 



1 1 — 

exact 

Hartree-Fock 




FIG. 20: Electron density for singlet states with different 7. 
Full line corresponds to exact results and dashed line to the 
Hartree-Fock approximation. Parameters as in Fig. 

In the case of triplet two-electron states, the single 
electron states are different, \<p) and \(p). Choosing these 
to be orthogonal, we get 



E- (i,i)t 



10) 



(CIO) 



The energy is now 



{il>\H\1>) = H#i|^<^> 

+ (^\H 1 \(p)( ( p\cp) + ^J2 U '' 



(Cll) 



and the norm is 



(ip\ip) = (<p\<p) (<p\<p) 



(C12) 



The system of equations for the coefficients ipi (and 
equivalent for (pi) is 

(0 \ciHi\ ifi) + ^ u ii l^il 2 W _ E Uijf>*f>i(pj = E h f(fii, 
3 i 

(C13) 



where 



E M = E + E. 



(l) _ (<p\Hi\<f) 



(C14) 



Eq. ( p!3 ) is a single-particle tight-binding Schodinger 
equation with Hamiltonian 



E 



i 

with potential 

h = €i +^2 \<pj\ 2 

3 

and renormalized hoping parameters 

iij = tSj^ + 2J Uijip*(pi. 
3 

The energy of the triplet bound state is then 



(C15) 



(C16) 



(C17) 



E = Ei 



hf ' 



hf ' 



El 



(i) 



(Ci8) 



In Fig. £l] the electron density for singlet (a) and for 
triplet (b) states are shown for different 7. Other param- 
eters are taken as in the case of the longer bulge, Fig. ||. 
As discussed in the text relating to Fig. ||, Hartree-Fock 
approximation for the singlet is less reliable since the 
Coulomb repulsion is stronger due to both electrons be- 
ing in the same state. Indeed, for 7 ~ 0.7 the Hartree- 
Fock approximation does not yield the bound state found 
in the exact result. 




-1.5 -1 -0.5 0.5 1 1.5 

z I a, 



1.5 -1 -0.5 0.5 1 1.5 
z/ a, 



FIG. 21: Electron singlet (a) and triplet (b) state density for 
various 7. Parameters are as in Fig. o. 
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